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ABSTRACT 

Important open questions in cosmology require a better understanding of the Big Bang 
singularity. In string and matrix theories, light-like analogues of cosmological singularities 
(singular plane wave backgrounds) turn out to be particularly tractable. We give a status 
report on the current understanding of such light-like Big Bang models, presenting both 
solved and open problems. 



1 Introduction 



The last two decades have witnessed enormous progress in our understanding of the com- 
position and evolution of the universe. One of the remaining challenges is to understand 
how the very early universe reached a nearly homogeneous, nearly flat state with a specific 
spectrum of density perturbations consistent with present observations. 

The most popular explanation is that the very early universe underwent a period of 
inflation [lj. If one assumes that inflation started and lasted long enough, it is able to 
explain the flatness and homogeneity of the universe. It also solves the monopole problem. 
The greatest success of inflation is that "simple" (single-field, slow-roll) inflationary models 
predict nearly scale-invariant, nearly Gaussian adiabatic density perturbations |2J. These 
are the seeds of large scale structure and are visible as temperature anisotropies in the 
cosmic microwave background. One may wonder, though, how the universe emerged in 
a state that allowed inflation to start. In other words, how was a suitably fine-tuned 
initial state selected? In particular, in general relativity, inflationary solutions are past 
geodesically incomplete (under a certain assumption that excludes a contracting phase in 
the past) [3]. The question should then be asked whether singularity resolution in a more 
fundamental theory puts constraints on which inflationary models are allowed. 

Alternatives to inflation include the cyclic universe [3]. The ekpyrotic mechanism (ultra- 
slow contraction) generates a spectrum of perturbations very similar to that of inflation, 
but in a contracting universe [5]. In general relativity, the transition from a contracting to 
an expanding (spatially flat) universe requires going through a singularity [6]. At present, 
it is unclear whether such a transition is possible and whether perturbations would go 
through essentially unchanged. The answer will have to come from a theory beyond general 
relativity. 

As we have argued, for the inflationary universe, and even more for alternative models, 
it is important to try and understand the big bang singularity. The work described in this 
review is motivated by several fundamental questions. Can we describe the big bang itself? 
How do space and time emerge from the big bang? Is it consistent to have a contracting 
universe before the big bang? Does the universe have a natural initial state, and if so, 
does it lead to inflation? String theory provides a short-distance modification of Einstein 
gravity, which is hoped to resolve space-time singularities. Existing formulations of string 
theory depend sensitively on the class of space-times one works with, and we will have to 
distinguish several classes of models. 

Perturbative string theory requires the background space-time to be specified from the 
onset. This background space-time has to satisfy supergravity equations of motion with an 
infinite number of corrections expanded in powers of a', the inverse of the string tension. 
In the high curvature regime, which necessarily accompanies singularities, all these a'- 
corrections generically become equally important, and the background equations of motion 
generically become completely intractable. One exception is provided by orbifolds, obtained 
from manifolds by discrete identifications [7] . Orbifolds contain new sectors of closed strings, 
namely "twisted" closed strings, which on the covering space connect a point and its image 
under a discrete identification. The rules of perturbative string theory on orbifolds are 
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inherited from those on the covering space. If the discrete identifications have fixed points, 
orbifolds can be singular. It is known that static orbifold singularities are resolved in per- 
turbative string theory, precisely thanks to the inclusion of twisted closed strings becoming 
light near the singular point. The hope 10 years ago was that time-dependent orbifolds (see 
[SI El QUI ED] for reviews) would lead to simple examples of cosmological singularities re- 
solved within perturbative string theory. It turned out, however, that at least the simplest 
models were plagued by divergences related to large backreaction and invalidating string 
perturbation theory (T2J [131 EH US] 13 It is worth noting that, in perturbative string theory, 
lightlike orbifold singularities are equally problematic as spacelike ones. 

One can therefore turn to non-perturbative formulations of quantum gravity, such as 
the AdS/CFT-correspondence [18] or matrix [19] or matrix string [201 I2U [22] theories. 
It should be said that those frameworks do not directly deal with the usual singularities 
discussed in cosmology. The AdS/CFT-correspondence requires the space-time to exhibit 
an AdS structure, at least asymptotically, while matrix and matrix string theories rely 
on the presence of a light-like isometry. None of these properties apply to the singular- 
ities naturally emerging in the context of classical cosmology: Friedmann, Kasner and 
Belinsky-Khalatnikov-Lifshitz space-times. Nevertheless, any progress on resolving light- 
like or space-like singularities, even in a not directly realistic model, would be very welcome, 
as it would point to mechanisms by which quantum gravity can in principle resolve cosmo- 
logical singularities. 

One clearly needs a compromise between the classes of space-times tractable within 
contemporary string theory and those relevant in cosmology. A number of directions of 
research have been proposed in this vein, and in this review we shall concentrate on one of 
them: analysis of light-like analogues of the usual singular cosmologies. For other models, 
we refer to the reviews [SI El QUI El [221 [21] • Within the class of light-like singularities, we 
will focus on studies in perturbative string theory and matrix (string) theory. For studies 
of light-like singularities using the AdS/CFT correspondence, see the review [25] . 

The prototype metric of the kind we are going to consider can be written as 



When fj, is diagonal and depends on x + through power-law functions, this looks very much 
like a Kasner solution, except that the dependences are on the light-cone time x + rather 
than on the usual time. When \x is proportional to the unit matrix with a power-law 
dependence on x + , the metric looks like a light-like version of flat Friedmann cosmologies. 
Furthermore, ([T]) actually arises when a Penrose limit is taken around a light-like geodesic 
hitting cosmological singularities [26] . 

The metric ([I]) is of the plane wave type, i.e., it describes a strong plane-fronted grav- 
itational wave (which is a non-linear generalization of the familiar linearized gravitational 
waves in Minkowski space-time). Such plane waves are most conveniently analyzed in the 

^ee, however, [5] for a suggested resummation of divergences by working in the eikonal approximation, 
[16] for a proposed resolution of an orbifold singularity in terms of orientifolds, and [17] for a closely related 
model in which the singularity is replaced by a phase with a condensed winding tachyon within perturbative 
string theory. 
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so-called Brinkmann coordinates, with their metric given by 



ds 
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2dx + dx~ + K ij {x + )x i x j (dx+) 2 + (dx { ) 2 , 



(2) 



where K,ij(x + ) represent the profiles of different polarization components of the wave. In 
pure gravity, Kij(x + ) is constrained by Ku = 0, giving the same number of polarizations 
as in linearized theory (a traceless symmetric tensor in D — 2 dimensions, with D being 
the number of dimensions of space-time). If a dilaton is present, Ku does not vanish and is 
related to the dilaton, which gives an additional independent polarization component. By 
a ^-dependent rescaling of x l , §Z§ can be brought to the so-called Rosen form, in which the 
metric only depends on u, making the planar nature of the wave front manifest. The Rosen 
form is implied in (Tj[|). (It is prone to coordinate singularities, however, and often avoided.) 

Plane wave space-times possess a number of remarkable properties in the context of 
quantum gravity theories. In perturbative string theory, they are known to satisfy back- 
ground consistency conditions to all orders in a' [27]. In other words, if supergravity 
equations of motion (zeroth order in a') are satisfied, all the higher order corrections will 
vanish automatically. This property comes from a special structure of the Riemann tensor 
in the plane wave background, and it allows perturbative string analysis of even highly 
curved plane waves. At the same time, the worldsheet theory of strings in plane wave 
backgrounds turns out to be especially simple (and linear, when the light cone gauge is 
imposed). Similarly, the light-like isometry needed for formulating matrix models is also 
present in plane wave geometries, so a matrix theory description of quantum effects in these 
space-times can be given. 

It is most natural to start by studying ([I]) with proportional to the unit matrix (i.e., 
an isotropic space). Such space-times can be thought of as light-like analogues of Friedmann 
cosmologies. By the equations of motion, such space-times cannot be empty: one needs to 
add matter to compensate for the curvature of the plane wave. A natural choice for this 
additional matter is the dilaton, a scalar always present in string theories and supergravities. 
If one starts with ([1]) and takes /ijj to be proportional to 5ij and depending on x + as a power 
law, the corresponding Brinkmann form metric ([2]) can be written as 



(where d = D — 2 is the number of i- indices). Thus, for negative k, the string coupling 
blows up at the singularity (x + = 0), invalidating string perturbation theory. One 
would then expect a perturbative approach to be of little use for addressing the question of 
singularity transition in that case. It can still be applied, however, if k > 0. (Singularities 
with positive k arise as Penrose limits of power-law Big Bang singularities, with the scale 
factor of the universe proportional to positive powers of times; in the same way, k < 
corresponds to Big Rip singularities.) 




(3) 



In turn, the equations of motion determine the dilaton to be [28] 

kd 

(f) = + C X + H In x + 



(4) 
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It is often thought that quantum-gravitational effects should naturally resolve singular- 
ities in some way. This would certainly be desirable, but is not so in our explicit examples. 
The singular plane wave ([3]) enters various quantum gravity constructions as the back- 
ground. The singularity is then always there at x + = 0, at least asymptotically, even 
though locally the geometry may be altered (or even dissolved by non-geometrical states). 
Mathematically, the singularity appears as explicit singular time dependence in the Hamil- 
tonians of string and matrix theories in the plane wave fl3]). One then has to understand 
how to deal with such singular time dependences. 

Free string propagation on ([3]) was studied in [2H]- In particular, it was suggested in 
that publication that the question of propagation across the 1/(X + ) 2 singularity in the 
metric can be addressed by employing analytic continuation in the complex X + -plane. 
Subsequently, in (29], another principle was proposed, which we motivate nextJl 

In string and matrix theories, it is necessary to satisfy stringent consistency conditions 
in order to maintain fmiteness and anomaly cancellation. In perturbative string theories, 
this question is very well studied, and it is known that the space-time background has to 
satisfy a'-corrected supergravity equations of motion in order for the theory to be well- 
defined (as already mentioned, the a' corrections are absent if one is working with plane 
waves, hence satisfying plain supergravity equations of motion is sufficient). For matrix 
theories, similar restrictions arise from considering ft-symmetry of the D-brane action [31 J, 
though the question does not appear to have been studied systematically. In any case, one 
would expect that the handling of the singularity should be subject to rigid constraints 
(given that even for smooth space-times one encounters rigid constraints). 

There is one approach to handling plane wave singularities that automatically takes 
benefit of what is known about smooth space-times and applies it to the singular limit. 
Namely, one can consider (j3J) as a limit of smooth metrics of the type (J2J), do the relevant 
computations, and take the singular limit at the end. Then, for any resolved space-time 
(EJ) consistency of string theory is guaranteed if (E} satisfies the supergravity equations of 
motion (without any further conditions). It is then natural to assume that the singular 
limit will likewise be a consistent string theory, provided that this limit exists. 

Even with these specifications, there are many ways to resolve ([3]). One class of resolu- 
tions appears to be very special however. The background (|3]) possesses a scaling symmetry 
and does not depend on any dimensionful parameters. It is natural to demand that this 
symmetry should be recovered when the resolution is removed. This will happen if the 
resolved plane wave profile does not depend on any dimensionful parameters other than the 
resolution parameter e. In this case, on dimensional grounds, 




(5) 



The limit ([3]) will be recovered if 




(6) 



2 Recent work on closely related models appeared in [50] . 
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for large values of 77, with some b > 2. 

In this review, we shall first concentrate on the perturbative string analysis of space- 
times 03]), and then discuss how this class of backgrounds can be treated in matrix and 
matrix string theories. 



2 Perturbative strings in singular plane wave backgrounds 
2.1 The light cone Hamiltonian and WKB solutions 

String worldsheet fermions are free in plane wave backgrounds [52] • We shall therefore 
concentrate on the bosonic part of the string action, given by 



J dt ^ do^-g (g ab G, u d a X»d b X» - ~atR®$) . (7) 



Ana 

We choose light-cone gauge X + = a'p + r and gauge-fix the worldsheet metric, 

det(g ab ) = -1, d a g au = 0. (8) 

After some algebra and in units a' — 1, one obtains the light-cone Hamiltonian as a sum 
of time- dependent harmonic oscillators: 

00 d 
n=0 1=1 

H K = + ±n(t/ £ )^a>!, (10) 



where X n are Fourier transforms in the cr-coordinate: 

X% a) = X l (t) + (cos {jut) X l n {t) + sin {jut) Kitj) . (12) 

ra>0 

The Hamiltonian (|9]) is quadratic and the solution to the corresponding Schrodinger 
equation can be found using WKB techniques, which are exact for quadratic Hamiltonians. 
The solution can be written as 

<P n {t;X) = A(ti,t)exp {iS cl , n [Xl n MKA) , (13) 

where S c ^ n [X\ n , ti\X^, t] is the "classical action" evaluated for the path going from X\ n at 
the time t\ to X % n at the time t, 

s cl [x\ n Mxit}= [\t>(^-(n* + ±n(*))^\. (id 
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If A n (ti, t) satisfies 

cP 1 
-2-An(ti,t) = Mti^g^S^Xi^ttlXlt], (15) 

then (fT3l) satisfies the original Schrodinger equation exactly. 

Up to normalization, a basis of such solutions, labelled by the initial condition X l n (ti) = 

x iw is s iven b y i 33 J 



II Jr-x *** \-^r £ KKnf^C - {Kfd t2 C + 2Xj >n x;] ) , (16) 



where C(ti,t 2 ) (suppressing the index n) is a solution to the "classical equation of motion" 
for the time- dependent harmonic oscillator Hamiltonian (11 II) : 

dlC(h,t 2 ) + (n 2 + \n(t 2 /e)) C(t u t 2 ) = 0, (17) 



with initial conditions specified as 

C(ti,t 2 )|t 1=t2 = 0, d te C(ti,t 2 )k =t2 = 1. (18) 

We shall refer to C(ti,t 2 ) as "compression factor". To derive the singular limit of the 
wavefunction ( ITBl it is sufficient to study the singular limit of ( fTTlfTBl . 



2.2 The singular limit for the center-of-mass motion 

For the n = mode, we obtain as the "classical equation of motion" 



X + \fl(t/e)X = 0. (19) 



f 2 



We need to study the e — > limit of the solution that obeys the initial conditions 

X(ti) = 0, X(tt) = 1, h < 0. (20) 

The singular limit of solutions to this equation has been analyzed in [33]. Performing 
a scale transformation Y(rj) = X(r]e), with rj = t/e, removes the e-dependence from the 
equation, leaving 

rj 2 

—Y + n( V )Y = 0. (21) 

This scale transformation is possible because our initial singular metric was scale-invariant 
and we have resolved it as in ([5]) without introducing any dimensionful parameters besides 
e. The existence of a singular limit is then translated [33] into constraints on the asymp- 
totic behavior of solutions to (12 ip . These "boundary conditions at infinity" are strongly 
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reminiscent of a Sturm-Liouville problem, and it is natural that a discrete spectrum for the 
overall normalization of Q is singled out by imposing the existence of a singular limit. 

For the specific asymptotics of our resolved profile ([6]), it can be shown [33J that, in the 
infinite past and infinite future, the solutions approach a linear combination of two powers 
(denoted below a and 1 — a, with a being a function of k, cf. (jSJIH]))- This power law behavior 
simply corresponds to the regime when the second term on the right hand side of ([6]) can 
be neglected compared to the first. It is then convenient to form two bases of solutions, 
one asymptotically approaching the two powers (dominant and subdominant) at r\ — > — oo, 



n-(??) = M a -+o(M a -), 

and another behaving similarly at 77 — > +00 



y 2 -(v) = \vr a - + o(\vr a -) 



Yx+(v) = \v\ a+ + o(M a+ ), Y 2+ ( V ) = M 1 -^ + od^n 



where a± is given by 



a± 



1 

2 + 



4 ± 



(22) 



(23) 



(24) 



(We are temporarily assuming that k can take two different values k± for the positive and 
negative time asymptotics, a possibility that will be discarded shortly.) The two bases are 
related by a linear transformation: 



Y2-(V) 



Q 



Y 1+ (v) 
Y 2+ (v) 



(25) 



where Q is a 2 x 2 matrix whose determinant is constrained by Wronskian conservation as 

W[Y^,Y^} = ^[F 1+ ,F 2+ ]detg. (26) 

The singular limit has been rigorously considered in [33J, but the results can be under- 
stood heuristically from the following argument [29J. Imagine one is trying to construct a 
solution Y to (|2~T1) satisfying some (e-independent) initial conditions at 771 = tx/e < 0. This 
solution can be expressed in terms of Yi_ and Y 2 - (a complete basis) as 



Y = C X Y X - + C 2 Y 2 . 



(27) 



Since the initial conditions are specified at rjx = tx/e, the asymptotic expansions (1221) are 
valid. There needs to be a non-trivial contribution from both Yi_ and Y 2 _ in the above 
formula in order to satisfy general initial conditions. Hence, the two terms on the right 
hand side should be of order 1. Therefore, we should have 



C 1= 0(e a -), C 2 = 0{t 



l-a_ 



(2? 



If we now apply ([25]) and (|23|) to evaluate Y at a large positive 77 = t 2 /e, the powers of e in 
C\ and C 2 will combine with the powers of e originating from Y" 1+ and Y 2+ and yield 



Y(t 2 /e) = Q n t a 2 + 0(e a - a +) + Q l2 t\- a+ 0{e a 



+Q2it 2 + 0(e 1 - a - a +) + Q 22 tl~ a+ 0(e a+ ~ a -). 



(29) 
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Since a + and a_ are greater than 1/2, this expression can only have an e — > limit if a + = a_ 
(i.e., k + = k_) and Q 2 \ = 0. The latter condition implies that the overall normalization 
of the plane wave profile will generically lie in a discrete spectrum, dependent on 

the specific way the singularity is resolved, i.e., the detailed shape of f2(r/). A particular 
exactly solvable example for this discrete spectrum (there called "light-like reflector plane" ) 
has been given in [31]. With Q21 = and det Q = —1, the matrix Q can be written as 



Q 



q q 
-11 



(30) 



with q being a real nonzero number (q does not affect the singular limit). For flat space- 
time we have q = 1 and for the "light-like reflector plane" of [31] we have q — — 1. In the 
singular limit, a basis of solutions is given by 



Y 1 (t) = (-ty, Y 2 (t) = (-t) 1 -", 

Y 1 (t)=qt a , Y 2 (t) = --t 1 ~ a , 

q 



t < 0, 
t > 0. 



(31) 



2.3 The singular limit for excited string modes 

Following our general discussion of free strings in plane wave backgrounds, the evolution of 
excited string modes is described by time-dependent harmonic oscillator equations 

^X(t) + (V + in(f/e)) X(t) = 0. (32) 

Solutions for the wavefunctions of the excited string modes can be expressed in terms of 
a particular solution to this equation C(ti,t 2 ) defined by (QZDCEBD- Hence, to analyze the 
singular (e — > 0) limit of the excited modes dynamics, it should suffice to analyze the 
singular limit of C{tx,t2). Because n 2 is finite, it is natural to expect that it does not affect 
the existence of the singular limit (governed by the singularity emerging from Q(t/e)). It 
can be proved that it is indeed the case for positive k [29]. 

The general strategy here is to analyze (132]) separately in the near-singular region (t 
close to 0) and the region where the e — y limit is regular. Away from t = 0, up to 
corrections vanishing as e is taken to 0, fl32|) can be approximated by 

X(t) + (n 2 + ^ X (t) = 0, (33) 

which is related to Bessel's equation. Around t = 0, one should expect that n 2 can be 
neglected, which leaves the equation for the zero-mode (already analyzed in the previous 
section) : 

X{t) + \{l(t/e)X(t) = 0. (34) 

More specifically, the separation into near-singular and regular regions should be organized 
as follows: 
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I II III 

u -t f t f 



We use t e to indicate a time that will approach zero in the singular limit as 

t e = e l ~ c t c , (35) 

with t staying finite in relation to the "moments of observation" t\ and t 2 . The number c 
(between and 1) should be chosen later as needed for our proof. 

One then can show that (for positive k) there is indeed a choice of c that makes deviations 
from the approximate equations (13311341) small in the appropriate regions. One can then 
construct the e — > limit of solutions to (132]) by taking approximate solutions satisfying 
fl33l) and (l34"j) . splicing them together and taking the e — ¥ limit in the end. Thus, one 
obtains exact expressions for the e — > limit of solutions to ( |32l . even though analytic 
solutions to ( 1321) at finite e cannot be given. 

We refer the reader to the original article |29j for detailed proofs, and here simply state 
the result: For k > 0, the singular limit of the excited mode evolution exists whenever it 
exists for the center-of-mass motion, and it is described by the following matching condi- 
tions: 



Y 1 (t) = V-tJ a -i/2(-nt), Y 2 (t) = V-tJi/2-a(-nt), t < 0, 

Y 1 (t) = qVtJar-i/2{nt), Y 2 (t) = -^-J 1/2 _ a (nt), t > 0, (36) 

where J u are Bessel functions. 

Note that there is a slight subtlety in the sense that convergence to the e — > limit is 
not uniform with respect to n (we have kept n fixed in our considerations). However, this 
does not affect the result as long as one is only interested in the limiting expressions at 
e = 0. More discussion is given in the original publication [29] . 

2.4 The singular limit for the entire string 

As we have seen in the previous section, for k > 0, consistent propagation of the string 
center-of-mass across the singularity guarantees that all excited string modes also propagate 
in a consistent fashion. This is not sufficient, however, to define a consistent evolution for 
the whole string, since even small excitations of higher string modes can sum up to yield an 
infinite total energy [27]. As we shall see below, the condition of finite total string energy 
(after the singularity crossing) turns out to be very restrictive. 

The total string excitation energy can be conveniently expressed in terms of the Bogoli- 
ubov coefficients for the higher string modes which can be extracted from (1361) as 

«n = ~- — —, r, 37 

2q sm(a!7r) 

.exp(-nra) + q 2 exp(ina) 
2q sin(a7rj 
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and they turn out to be independent of n. Here, a = y/1 — Ak/2. The total mass of the 
string after crossing the singularity is given by [27] 



Since the j3 n are n-independent, M can only be finite if /3 n — for all n. For k > 0, this 
cannot be achieved, since < a < 1/2 and q is real. (For k = 0, which is the case of the 
"lightlike reflector plane" of [31], all f3 n will vanish if q 2 = 1, which is satisfied automatically 
for any reflection-symmetric Q.) The implication is then that, if the singularity is resolved 
without introducing any new dimensionful scales, the singularity transition cannot be de- 
fined. (The alternative is having dimensionful parameters buried strictly at the singular 
locus, and an explanation of the possible physical origin of such parameters would be in 
order.) 

As mentioned before, in the case k < the string coupling blows up near the singularity, 
so that string perturbation theory is certainly not valid. Our considerations can be seen as 
a motivation to study these backgrounds in the context of non-perturbative formulations 
of string theory, to which we now turn. 

3 Matrix big bang models 

3.1 Time-dependent matrix and matrix string theories 

Matrix theory [Tj| is a non-perturbative formulation of M-theory in 11-dimensional Minkow- 
ski space-time. One way to derive it is by Discrete Light-Cone Quantization (DLCQ) [35J. 
A sector with N units of lightcone momentum is described by a quantum mechanics of 
N x N matrices, namely the dimensional reduction of (9+l)-dimensional SU(N) super- 
Yang-Mills theory to 0+1 dimensions. In matrix theory, (lightcone) time is built in, but 
space is an emergent concept, arising from a "moduli space" of flat directions. For this 
"moduli space" , and therefore space-time, to emerge, supersymmetry plays an essential role: 
without supersymmetry, quantum corrections would lift the flat directions. Compactifying 
M-theory on a circle leads to type IIA string theory. Matrix string theory [201 [22] is a 
non-perturbative formulation of type IIA string theory in lOd Minkowski space-time. The 
aspects of matrix and matrix string theory needed for our purposes are reviewed in detail 
in [TO]. 

The construction of matrix theories essentially relies on compactifying a light-like di- 
rection (which is a pre-requisite for discrete light-cone quantization). Due to discreteness, 
positivity and conservation of the light-cone momenta in such a compactified space-time, 
Hilbert spaces of quantum theories living in this space-time split into independent sectors 
(labelled by the value of the light-cone momentum), each of which correponds to a quan- 
tum system with a finite number of degrees of freedom. For matrix theories, this quantum 
system is given explicitly by finite N matrix Lagrangians. This is a remarkable simplifica- 
tion (which automatically renders the theory UV-finite, since there are no divergences in 
quantum mechanics). The drawback is equally grave, however, as reconstructing quantities 




(39) 



71 
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in an infinite space-time from the compactified version is highly non-trivial (it is believed 
to be possible, however, since for a sufficiently large compactification radius, an arbitrarily 
large laboratory can fit into the space-time with a light-like compactification). 

Compactifying a light-like direction is only possible in space-time backgrounds with a 
light-like isometry. This obviously excludes the case of ordinary cosmological space-times, 
as well as almost all familiar non-trivial solutions to the classical gravitational equations of 
motion. Strong plane waves of the type (pfl E]) are an exception though, which highlights once 
again their special status in quantum gravity and makes them an interesting laboratory for 
applying matrix theory methods to study time dependence and strong gravitational effects. 

One must keep in mind that reconstructing quantities in a decompactified space-time 
from their DLCQ analogues (explicitly given by the finite N matrix theories) may become 
more subtle in a time-dependent setting H The point is that, at best, DLCQ can provide 
information on quantities in a decompactified space that are measurable in a finite-size 
laboratory. This is not likely to pose a problem for the case of mild time dependences. 
But for the opposite extreme, namely space-times including singularities (which we shall 
comment on in the next section), linear dimensions of physical systems in an infinite space 
may blow up indefinitely. In that case, finite box dynamics (and hence DLCQ) will not 
be adequate to describe the evolution in a decompactified space. Whether or not these 
subtleties do arise has to be decided on the basis of careful dynamical considerations, which 
have not been carried out yet. 

We shall now give a summary of a few different versions of matrix and matrix string 
theories in plane wave backgrounds. For the original time- dependent Matrix Big Bang 
matrix string theory of [36J, the 10-dimensional geometry is asymptotic to the linear dilaton 




To construct the matrix string theory for the background fj4*0l) . one first lifts the background 
( I40p to 11 dimensions via the usual conjecture of type IIA/M-theory correspondence. The 
resulting 11-dimensional space-time is 



where y is a coordinate along the M-theory circle. This is followed by the DLCQ compact- 
ification of the light-like f-coordinate, interpreted as the M-theory circle of an "auxiliary" 
type IIA string theory. A T-duality [IB] then relates the resulting theory of DO-branes on 
a compact dimension, i.e., a BFSS-like matrix theory with a compactified dimension, to a 
more manageable theory of wrapped Dl-branes. This procedure has been carried out (in a 
slightly different but equivalent way) in [36] and has been reviewed in [10] and [27]. The 

3 We thank David Kutasov for drawing our attention to this issue. 

4 The original set-up of [35] has been later extended in various directions [371 EH EH1 SOI HH S21 SSI SH 
I45[ 1461 [47] ; in particular, a systematic generalization of the analysis to more general singular homogeneous 
plane- wave space-time backgrounds has appeared in |47j . 



configuration^ 





(40) 



= -Qy 



ds 2 = e 2Qy+ / 3 (-2dy + dy- + (dy 1 ) 2 ) + e~ iQy+ '\dy) 2 , 



(41) 
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resulting matrix string action is 

s = {^{d,x^ + e T ij>e + ^Fl- g l h Ax\x^f + g YM e T 1 Ax\e}^ , 

(42) 

with the Yang-Mills coupling gyM related to the worldsheet values of the dilaton: 

e -0(?/+(r)) e Qr 

9ym = k—j = 7T1 — • ( 43 ) 

2irl s g s 2irl s g s 

A generalization of this set-up has been proposed [47]|. One can start with a 10- 
dimensional power-law plane wave: 



i\2 



ds 2 st = -2dy + dy- + g^dyW = -2dy + dy- + J2(y + ) 2m *(dy 

i 

= -2dz + dz- + ma ^~ l \ z a f{dz + f + ^(^ a ) 2 , ( 44 ) 

e 2< ^ = (y+) 3fe /( b + 1 ) = ( z + ) 3b/(b+1) . 

Here, the first and the second line represent the Rosen and the Brinkmann form of the same 
plane wave, respectively. In order for the supergravity equations of motion to be satisfied, 
one needs to impose [4T] 

J^K-lH-^y. (45) 

i 

The original background of j3S] can be seen as a b — > — 1 limit of the above space-time [47]. 
The 11-dimensional space-time corresponding to (1441) is 

ds 2 u = -2dudv + J2 u2n '( d y i f + u 2 \dyf 

(46) 



-2dudw + Y na(n " 1} (x a ) 2 (^) 2 + MLJl x 2 (du)2 + yr dx ay + , d 



x) 2 , 



with rii related to m; by 2mj = (2rii + b)/(b + 1). The usual formulation leads to a matrix 
string action, whose bosonic part is given, in the Rosen coordinates of (1441) . by 

S RC = f drdcr Tr ^—^gy^r]^^ 6 FupF^s — ^p a ^ g%i{j)D a X l D ^X^ 

+ \g YM g lk {r)g 3l {r)[X\X^[X\X 1 ]^ . (47) 

with the transverse metric g^ given by the first line of (144|) and the Yang-Mills coupling by 

e -0(y + (r)) r -36/2(6+i) 

9ym = — tti = • ^48) 

2nl s g s 2nl s g s 
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One can further transform this action to the Brinkmann coordinates of the original plane 
wave, given by the second line of (jSJ), to obtain [4"T] : 

Sbc = j drdaTr (-~g? M F? ff - \ (D T Z a D T Z a - D a Z a D a Z a ) 

+\g 2 YM [Z a , Z b ][Z a , Z h ] + l -A ab {r)Z a Z b ^ , (49) 

where A ab = diag{m a (m a — l)}/r 2 . The latter form of the action only differs from a SYM 
gauge theory with a time-dependent coupling by the term involving A ab . 

In [37], 11-dimensional (quantum- mechanical) matrix theories were introduced as sim- 
pler analogues of the matrix string theories of [36J. The relevant 11-dimensional (M-theory) 
background has the form 

ds 2 = e 2ax+ (-2dx + dx~ + (dx') 2 ) + e Wx+ (dx 11 ) 2 , (50) 

or, in terms of the light-like geodesic affine parameter r = e 2ax+ /2a, 

ds 2 = -2drdx~ + 2ar{dx i ) 2 + {2aT) l3/a {dx 11 ) 2 . (51) 

This metric satisfies the 11-dimensional supergravity equations of motion if the constants a 
and (3 are related as (3 = —2a, or (3 = 4a. The fact that these relations need to be imposed 
will not be relevant for what follows (it is essential, however, for the general consistency 
of the corresponding matrix theories). Since translations in x~ form an isometry of the 
above background, the usual DLCQ argument (proposed in [49j and adapted to the time- 
dependent case in |36j) can be applied. The result [37J is a matrix theory that can be 
expected to describe non-perturbative quantum gravity in space-times asymptotic to (|5Ti) . 
The bosonic and fermionic parts of the matrix theory action, respectively, have the following 
form: 

Sb = j rfrTr {^(D T Xr + ^^(DrX 11 ) 2 - ^{2arf[X\X^ 

_|(2« r ) 1+ ^[X i ,X 11 ] 2 |, ( 52 ) 
S F = Jdr {i6 T D T 6 - RVfaF& r >y i [X i ,e] - R{2arf' 2a e T ln [X 11 , 6}} . 
3.2 Singularity transition 

The models we have just formulated are aimed to describe quantum-gravitational effects in 
singular plane waves In Rosen coordinates (first line of (14"?)]) ). these space-times can 

be seen as a light-like analogue of Friedmann or Kasner cosmologies. It is then interesting 
to investigate what the theory has to say about the question of singularity transition. 

The general structure that emerges from our considerations is matrix (string) Hamilto- 
nians with singular time-dependences. This can be explicitly seen in (}4"9"]) as A ab is singular 
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at t = 0. Likewise, the factors of r in front of the time derivatives in fl52l) will become 
inverted when the velocities are replaced by momenta, and the corresponding Hamiltonian 
will contain an isolated singularity in its time dependence at r = 0. These isolated sin- 
gularities in the explicit time dependence of the matrix (string) Hamiltonians are directly 
inherited from the singularity in the background space-time. (Such singularities are also 
present in Hamiltonians (llOHlip in the perturbative string theory setting of the previous 
section, if one looks at the e — > singular limit directly without performing the singularity 
resolution.) 

One might have thought that quantum gravity should resolve singularities in some way 
and give a dynamical prescription for singularity transitions. Unfortunately, it is not so in 
our present setting. The underlying space-time singularities appear as explicit singularities 
in time dependences of matrix (string) Hamiltonians, and the evolution is not well-defined 
without additional prescriptions. 

Some general properties of Hamiltonians with singular time-dependences have been 
studied in [5U1 EI]. An important thing to understand is that there is a tremendous ambi- 
guity associated with defining the singularity transition in this context. Indeed, if one only 
assumes that the Schrodinger equation for the Hamiltonian with a singular time depen- 
dence is satisfied away from singularity, and finds a way to define the singularity transition, 
one can immediately create another singularity transition prescription with the same prop- 
erties. Namely, one can change the wave vectors by an arbitrary unitary transformation 
the moment they pass the singularity. The Schrodinger equation will still be satisfied away 
from the singularity, as it was before. One is, therefore, in a need of a physical (or at least 
heuristic) motivation for choosing a particular prescription for singularity transition. 

We believe that a very natural class of singularity transitions arises from considering 
geometrical resolutions of the singular plane wave (T4T)|) . Due to the functional arbitrariness 
of the plane wave profile, it is very easy to replace the singular profile in (T4T)j) by a resolved 
one in a way similar to ([5J16]). Maintaining the background space-time as a solution to the 
equations of motion is also essential from the standpoint of background consistency. En- 
forcing supergravity equations for the background appears to be related to the re-symmetry 
of the D-brane action [21] necessary for the standard formulation of matrix theories. Unlike 
the case of perturbative strings, the action is non-linear, so the singular limit may be much 
more difficult to investigate, and it remains an important open problem for the future. 

(Note the similarities of this setting to the AdS light-like cosmologies reviewed in [25J. 
There, one ends up with a gauge theory featuring singular light-like time dependences. The 
only difference in our case is that the singular dependences are on ordinary time. Similar 
non-linearities characteristic of gauge theories are present in both cases. In [22], a special 
class of plane waves is chosen for which the analysis simplifies due to conformal flatness. 
General p-brane-plane-wave solutions presented in [51] suggest a natural generalization of 
AdS light-like cosmologies featuring arbitrary plane wave profiles. In that setting, the 
problem of singularity transition is even more similar to what one encounters in Matrix Big 
Bang models.) 
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3.3 Near-classical late time dynamics 

Even though, in the context of time-dependent matrix and matrix string theories, novel 
physics is expected to emerge in the high-curvature regions of space time, it is also im- 
portant to understand how the near-classical space-time emerges away from the singularity 
when the curvature becomes small. The issue may be technically somewhat involved, since 
geometrical notions appear rather indirectly in the matrix formalism. Heuristic results ap- 
peared in [361 E2]; a systematic analysis has recently been carried out in [53], the results of 
which we now review Jj 

The problems of late time (near-classical) dynamics in matrix models can be understood 
by inspecting the action (I52"j) . The space-time backgrounds implicit in (152]) are supposed 
to feature a light-like singularity at r = and become progressively more classical at 
large r. Yet, the explicit time dependences in (152]) superficially become more steep, if 
anything, at large r. Additionally, there are no supersymmetries explicit in ([52]). Since 
supersymmetries are crucial for the free propagation of well-separated gravitons (and hence, 
a robust geometrical interpretation) in the flat space matrix theory, one should attempt to 
find an analogue of supersymmetry in (152]) that would enforce a similar type of dynamics. 

To address these important issues, we first note that the metric (1511) describes a plane 
wave and, with the coordinate transformation 

u = r, z l = V2^x\ z 11 = {2arfl 2a x l \ v = x~ + ± , (53) 

4ar 

it can be brought to the Brinkmann form 

ds 2 = -2du dv — 7- — .„ ,y — -du 2 + idz 1 ) 2 + (dz 11 ) . (54) 

(2au) 2 

This parametrization forces the metric to manifestly approach Minkowski space-time for 
large values of the light-cone time, which strongly suggests that the large time dynamics 
of the corresponding matrix theory will likewise approach the flat space matrix theory, if 
treated in appropriate variables. (As we shall see below, the convergence towards this limit 
is somewhat subtle, but the naive expectation will prove well-grounded.) 
The matrix theory corresponding to (1541) is given by 

S B = j drTr ji \{D T Z i ) 2 + (D T Z n ) 2 } - | [[Z\ Z>} 2 + 2[Z\ Z n ] 2 } 



a 



2l V) 2 -{fi 2 -2afi){Z 1X f \ (55) 



(2ar) 2 

S F = jdr {iO T D T 6 - R6 T ll [Z\ 6} - RQ T lxx \Z xx , 9}} . 

The action ( 155]) only differs from the flat space matrix theory by a term decaying as 
1/r 2 , thus one may expect that the late time dynamics will be approximated by the flat 



5 Another attempt to study late time (low background curvature) dynamics of the time-dependent matrix 
theories was undertaken in 1541 . 
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space matrix theory and admit the usual space-time interpretation. However, the decay is 
quite slow and one might be worried about whether it is sufficient to ensure convergence. 

To illustrate these worries, one may look at the straightforward example of a harmonic 
oscillator whose frequency depends on time as 1/t 2 : 



x + —x = 0. 
t 2 



(56) 



The two independent solutions to this equation can be given as t a and t 1_a , where a is a 
fc-dependent number. These two solutions are obviously quite different from a free particle 
trajectory, even though the equation of motion approaches that of a free particle at late 
times. The reason for this discrepancy is the slow rate of decay of the second term in (1561) . 

However, in a physical setting, one is only able to perform finite time experiments. That 
is, one has to specify the initial values x(t ) = x , x(t ) = v and examine the corresponding 
solution between to and to + T. The solution is given by 



x(t) 



x {l - a) - v t ( t 



I -2a 



to 



+ 



Voto — XqCL ( t 



1 - 2a 



to 



1-a 



(57) 



One can then see that x(to + T) = xq + v$T + O(T/t ), i.e., it is approximated by a free 
motion arbitrarily well if the experiment starts sufficiently late. 

It may be legitimately expected that the finite time behavior of the full time-dependent 
matrix theory given by ( T55|) will be approximated arbitrarily well by the flat space matrix 
theory at late times, just as in the above harmonic oscillator example. We shall now 
prove it by constructing an elementary bound on dynamical deviations due to a small 
time-dependent term in the Schrodinger equation. 

We start with the following Schrodinger equation: 



i^|$> = (JT + /(f)ffi)|*), 



(58) 



where Ho and H% are time- independent, and rewrite it in the interaction picture (with 
respect to H ): 



1$) = e - lHo{t - to) |o 



We then proceed to consider 



f( t yiH {t-t ) Hie -iH (t-t ) |^ _ 



(59) 



d_ 

dt 



j t mto) \m)+cc.) 



(60) 



\m) - mo)) 

= ~if(t) {mo^e^-^H^^ \£(t)) ~ c.c.) . 
Integrating this expression between to and to + T and making use of standard inequalities 
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for absolute values and scalar products, we obtain: 



t +T 

left, + T)> - |eC*d)> f = ~i J <&/(*) ((e^'^H.e-^-^ato) |£(f)> - cc.) 

to 

to+T 

<2 | ^|/(t)i V /(|e^o(*-*o) jf f ie -^a(*-*o)^( fo )» 2 v /(| e (t)»2 (61) 

to 



<2(max [t0it0+T] |/(t)i) / dtJ(ato)\e^Hfe-^\ato)} 



o 

Now, assume that f(t) approaches at large times and consider a fixed |£(^o)) = |£o) ( so 
we consider the evolution with fixed duration T of the same initial state |£o) starting at 
different initial times to). In this case, the first factor in the last line becomes arbitrarily 
small for large to, whereas the second factor does not depend on t . We then conclude that, 
for sufficiently late times, the finite time evolution of the state vector will be approximated 
arbitrarily well by = const, i.e., by the evolution with /(£) set identically to 0. 

It is then a simple corollary of the above bound that the time-dependent matrix theory 
dynamics becomes approximated arbitrarily well at late times by the flat space matrix 
theory, and, in particular, the supersymmetry is asymptotically restored (with all the usual 
consequences, such as protection of the flat directions of the commutator potential, and 
free graviton propagation). 

Similar considerations can be given for the matrix string case, though they are more 
involved and rely on derivations in the style of quantum adiabatic theory. We refer the 
reader to the original publication [53J. 



4 Conclusions 



We have reviewed some recent considerations of light-like singularities in string and matrix 
theories. In the presence of such singularities, these theories need to be supplemented with 
well-motivated prescriptions for singularity transition. A natural class of such prescriptions 
emerges from geometrical resolutions of the light-like singularities for which resolved space- 
times satisfy supergravity equations of motion. Important open questions include analysis 
of the singular limit for the full non-linear evolution of matrix theories. 
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